CERN-PH-TH /201 1-282 
NORDITA-201 1-101 



Inhomogeneous non-Gaussianity 

Christian T. Byrnes," Sami Nurmi,'' Gianmassimo Tasinato,'^ David Wands'^ 

"CERN, PH-TH Division, CH-1211, Geneve 23, Switzerland 
''NORDITA, SE-106 91, Stockholm, Sweden 
^Institute of Cosmology & Gravitation, University of Portsmouth, Dennis Sciama Building, 

Portsmouth, POl 3FX, United Kingdom 

■ We propose a method to probe higher-order correlators of the primordial density field through 
I the inhomogeneity of local non-Gaussian parameters, such as /nl, measured within smaller patches 
, of the sky. Correlators between n-point functions measured in one patch of the sky and fc-point 

functions measured in another patch depend upon the (n -I- fc)-point functions over the entire sky. 
^ ' The inhomogeneity of non-Gaussian parameters may be a feasiljle way to detect or constrain higher- 

order correlators in local models of non-Gaussianity, as well as to distinguish between single and 

■ multiple-source scenarios for generating the primordial density perturbation, and more generally to 
^SJ ' probe the details of inflationary physics. 



X: 



I. INTRODUCTION 



o 

. Primordial non-Gaussianity is a window onto the early universe, which can be used to distinguish between 

Qj different models for the origin of cosmological fluctuations [1-5]. For this reason, it is important to characterize 
as completely as possible the properties of observable quantities in a non-Gaussian distribution. In this paper, 
$H ' we discuss under what conditions local non-Gaussian observables can be inhomogeneous to a measurable level. 
0^ I Consider subdividing the observable sky of size into a large number of small patches of size £ <^ H~^, 

Ci • and measuring cosmological observables on each of these patches. The effective values of the cosmological 
parameters can differ from one patch to the other, since each small part of the sky experiences a different local 
. background. In particular primordial curvature fluctuations within a small patch are defined with respect to a 
^ ■ local background which can be decomposed into a cosmological background solution of the homogeneous field 
equations, plus the cumulative effect of random modes with wavelengths larger than the size £ of the patch 
■ (but smaller than the size H^^ of the observable sky). Since different patches of our observable universe are 
, affected by different contributions from the random long wavelength modes, each patch is defined on a different 
background configuration. As a consequence, quantities characterizing the properties of n-point functions of 
curvature fluctuations can be different in one patch or the other. Indeed the effect of long-wavelength modes 
on the small-scale gravitational collapse of dark matter halos is described by the peak-background split in 
astrophysics and is used to explain the biased clustering of dark matter halos [6]. In a non-Gaussian field the 
long-wavelength modes also modulate the amplitude of the short-wavelength density perturbations which can 
lead to a distinctive scale-dependent bias which is a powerful probe of primordial non-Gaussianity [7]. Of course 
there is also a cosmic variance between patches, even if there are no non-zero higher order correlators. However 
this effect will be small as long as there are a large number of independent measurements in every patch. 
] Correlators between n-point functions measured in one patch, and fc-point functions measured in another 

patch, depend on (ri + fc)-point functions in the entire sky. As an example, we will show that the autocorrelation 
of /nl between different patches depends on a class of parameters controlling 5- and 6-point functions, that 
in the single-source case reduce to squares of trispectrum parameters. This implies that inhomogeneities in 
quantities characterizing non-Gaussian observables can be measurable, since at present we have only weak or 
even no constraints on parameters controlling higher-order point functions in the full sky. Moreover, although 
the optimal way to constrain higher n-point functions in the full sky would be through a direct analysis of that 
correlator, in practice this is exceedingly time consuming and computer intensive. This problem will become 
significantly more pressing when the Planck satellite data becomes avalailable in early 2013 [8], since it has a 
much higher resolution than the WMAP satellite [9] . Only during the last few years have the first constraints 
on the trispectrum been made [10-14], and there are not yet any constraints on higher n-point functions. 
Therefore, our method of relating n-point functions to inhomogeneities of lower-order functions should also 
provide a feasible and practical way of constraining higher-point functions in the full sky, even though the 
method is suboptimal. There is a clear theoretical interest in constraining higher-order correlators: they may 
be much larger than the correlators which have already been constrained, and their detection would provide 
important additional information for characterizing the origin of primordial fluctuations. Refs. [15-17] consider 
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the amplitude of correlators of higher order than the trispectrum for local models of non-Gaussianity. See [18] 
for an analysis of non-Gaussianity not characterized by a local shape, that investigates 10-point functions and 
beyond. 

The inhomogeneities of non-Gaussian quantities that we analyze in this paper are distinct from the anisotropics 
of inflationary observables, arising when background vector fields (perhaps of curvaton nature) are turned on 
[19 27]. We are going to consider primordial non-Gaussianity due to scalar fields that do not induce statistical 
anisotropy. Conversely, generic models of inflation with vector flclds may produce statistical anisotropy but 
do not induce inhomogeneities of inflationary parameters. The inhomogeneities that we consider here are also 
diff'erent from the scale-dependence of the non-linearity parameters, that describe a variation with size rather 
than with position of patches [28]. 

The plan of this paper is as follows: We first focus on single-source scenarios, in which only one field generates 
the curvature perturbation. In Sec. II A we consider correlators of ^ and the power spectrum, in Sec. II B of /nl 
and in Sec. II C possible tests of single versus multi-source scenarios. In Sec. Ill we focus on the more general 
formula of multi-source models, before going on to analyze in Sec. IV when the inhomogeneities can be large 
enough to be observable. Finally we conclude in Sec. V. 



II. STATISTICAL INHOMOGENEITIES IN SINGLE-SOURCE MODELS 



Let us start by considering single-source models where the curvature perturbation ( arises from fluctuations 
of a single scalar fleld. The single-source scenarios that we have in mind do not necessarily correspond to 
single-field models of infiation. We have in mind set-ups in which more than one scalar field may be present 
in the system, but only one of them, which we dub a, is responsible for generating the primordial curvature 
perturbation. Examples are the curvaton model [29, 30], or the modulated reheating scenario [31, 32]. These 
models are particularly interesting since they are capable of generating large non-Gaussianity of the local type, 
potentially observable by the Planck satellite. 

The statistics of the curvature perturbation in general depend on the size of the patch we consider. This is 
due to long- wavelength fiuctuations which contribute to local background quantities in any patch smaller than 
the horizon. The curvature perturbation within a patch of size £ is given by [33-37] 

0(x) = N'{at)Sta{^) + iiV"(a,)5,a(x)2 + ^N'" {ae)Sea{^f + . . . , (1) 

where (T£ denotes the background field value in the patch. It consists of the classical homogenous solution d-{t), 
and of fiuctuations da^ with wavelengths greater than the patch size, k < a(r^ (here a denotes the scale factor). 
Nia) = J H dt denotes the number of e-foldings from an initial time ti, soon after horizon exit during inflation 
of the modes of interest, to some final time (for example, during the primordial radiation-dominated era) at 
which ( has frozen to its final, constant value. The fluctuations 5e<T{x) include modes ranging from k = a£~^ 
to ki = QiHi, the latter corresponding to modes that exit the horizon at the initial time ti ^ : 

^^"^")= / (^^^''■'^"''' • 

a^-i<fe<fei 

We use a top-hat window function to select the modes within this momentum interval. Different choices of the 

window fimction should not significantly alter our results. (5£(t(x) consists of superhorizon modes at the initial 
time, (JcTk, which can subsequently be treated as a classical, Gaussian random field with the two-point function 
given by ^ 

{Sa^Sak'} = (27r)35(k + k')^ • (3) 

Here the brackets denote ensemble averages. 

We are interested in computing ( in patches smaller than the observable universe, i <^ H~^, see Fig 1. The 
local background field value in a patch centered at some position x is given by 

/dk' / 
J^:^^'^ '''^^k' = O-Jf-i + A^CTx , (4) 

aH<k<ae-^ 



Notice, however, that this cut-off docs not lead to an additional de])endciice on ti of the quantity Sicr{x), besides the one already 
contained in 5(7^.. Indeed, the upper limit ki is meant to characterize a UV cut-off. Pushing this cut-off to a slightly larger fixed 
scale, k-av > ki would remove the additional dependence on ti, without qualitatively changing our results. This since subhorizon 
modes are expected to provide only subdominant contributions to the evolution on superhorizon modes. 
^ We assume that the scalcir field cr has canonical kinetic terms, and is chajracterized by slow-roll dynamics during inflation. 
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where cr^f-i denotes the background field in the entire observable universe. A^CTx is comprised of fluctuations 
(5(Tk with wavelengths greater than the patch size, k < a£~^, which do not average out when computing spatial 
averages over the patch i. For a patch i located at a position x, Aga^ therefore acts as a constant, local 
background. We emphasize that we use the label x to indicate quantities evaluated in a patch at position x, 
and we consider them not as functions of the coordinate x. Denoting the spatial average over the patch by 
(. . .)e, we in general obtain the non-vanishing result 

Aea^ = {Aia^)i ^ , (5) 

which depends on the location x of the patch. On the other hand, when £ <C H^^, the spatial average of A^ctx, 
when computed in the entire observable universe H~^, vanishes: (A£crx)ij-i = 0. This is because A^Cx in our 
approximation depends linearly on the fluctuations dai^, which have a vanishing spatial average over the full 
sky. By the ergodic theorem, and provided that the quantity i ■ H is suSiciently small, the average of A^cTx 
computed over the full sky coincides with the ensemble average, (A^crx)/f-i = (A^CTx) = 0. 




FIG. 1: Schematic diagram to explain the measurements we are proposing. The large shaded area is the full sky (Hubble 
volume) while the two small regions of size I <^ are two examples of the small patches in which observations are 
made and then correlated. 

The consequences of long- wavelength fluctuations have been extensively studied in the context of the infrared 
growth of inflationary correlators. It has been shown that, for sufficiently local observations, a suitable shift of 
background quantities is able to remove infrared divergences associated with adiabatic field fluctuations [38- 
41], see also [42-44]. On the other hand, infrared-enhanced effects associated with long wavelength isocurvature 
fluctuations cannot be removed by shifts in background quantities, and have the opportunity to provide sizable 
contributions to inflationary obscrvablcs [45]. In our analysis, we will consider correlators between cosmological 
observables measured on distinct subhorizon patches within the observable universe, in set-ups where long- 
wavelength modes generate sizable statistical inhomogeneities. This effect was pointed out already in [46], 
but for the specific case of single-field slow-roll infiation where inhomogeneities arc slow-roll suppressed. Our 
treatment applies to generic single and multiple source models, and we find that the inhomogeneities can become 
large in scenarios characterized by observable non-Gaussianity. We will only consider the inhomogeneities arising 
from scalar perturbations. Tensor modes, included in the analysis of [38, 39], should only generate subleading, 
slow roll suppressed corrections to our results which we neglect. 

After this general discussion aimed to define the set-up, we move on discussing in more detail small scale 
statistical inhomogeneities of quantities characterizing the properties of the curvature perturbation. 



A. Correlators of one and two-point functions 

The two-point function of the curvature perturbation (1) in a patch of size i <C is given by 

(0(yi)0(y2))£ = N'{ae,^f{Seaiy,)5ea{y2))i + 0(/SlC*) (6) 

If we re-express this with respect to the background field in the entire observable universe, ajj-i, using Eq. (4) 
this becomes 

{Q{yMy2))e = N'{aH-^f{Sea{yr)Sea{y2))e (l + '^^'^^^^^^ ^^^^ + ^(/nlC') , (7) 

where x denotes the location of the patch and yi , y2 are coordinates inside the patch. Recall that brackets ()^ 
denote spatial averages computed over the region I. The long-wavelength contribution A^ctx breaks translational 



4 



invariance: the two-point function depends not only on the separation of the points, |yi — y2|, but also on the 
location of the patch x. The long-wavelength modes therefore generate statistical inhomogeneity. For the case 
of 2-pt functions, this fact was also pointed out in [46] and, in more generality, in [48]. 

Taking a Fourier transformation of (7) within the region ^, i.e., integrating over the y-coordinates and treating 
X as a constant, we find a relation between the spectra of Q and Ch-i 

A,n(x)=n.-n.-i = ^Vc^_,hi.{aH-^)N'{aH-^)Aea^ + 0{^„IN'vl'^) + 0{fl,^Vl) , (8) 



where the full sky power spectrum is 



n.-.=^WO^(f) , (9) 



and /nl describes the amplitude of the primordial bispectrum relative to the square of the power spectrum, 
which is given in the JA/'-approach by [37] 

5 N" 

/NL=g-^- (10) 

A£7-'^(x) measures deviations of the local spectrum, measured in a patch of size I, from the global power 
spectrum characterizing 2-pt functions of perturbations over the entire observable universe. 

Since A^Pf (x) is proportional to A^ctx, it is a Gaussian field with zero mean over the full sky for t <C 



with two-point function 



(A^Pf(xi)A^Pf(x2))H-i = V^^_^ (y/NL(aH-i)^ iV'(aH-i)'(A^a(xi)A^a(x2))H-i (11) 



where denotes the spatial average computed over the entire observable universe, which we assume 

coincides with the ensemble average (...). Hence we will drop the label from the angle brackets in what 
follows. We denote Aa; = jxi — X2I, ke = ai~^ and fchor = aH. The integral can be evaluated as 



12 ^^ / ,j , / sin(fc'Aa;) 



F{Ax,khor,ke) = dk 



, sin(fc'Aa;) 



"'hor 

... . . .,, . . smjkeAx) sin(A;horAa:) 

= ci(fc,Ax) - ci(fc,o.Ax-) - + fe^^^^^ , (12) 

with ci(.x) the cosine integral function. Notice that, as expected, the correlators between two-point functions 
evaluated in different patches of the universe are proportional to (|/nl) = tnl, which in single-source models 
is a parameter characterizing the four-point function measured in the entire sky.^ In the limit Aa; 0, the 
result takes a particularly simple form 

F{Ax = 0, fchor, ke) = In = -In {iH) , (13) 

corresponding to the well-known logarithm associated with IR-enhancements. For kj^ <^ Ax -C k^^^ it can be 
approximated by 

^^(Ax, fchor, fc£) « 1 - 7E - In(fchorAx) + O ((fc^Ax)-!) (14) 

where 7e ~ 0.58 is the Euler constant. A comparison between the exact result (12) and the estimate (14) is 
shown in Figure 2. 

The magnitude of the statistical inhomogencitics is described by the variance of AfVc^. Using COBE normal- 
ization Vc = 2.4 X 10"^ and the bound |/nl| < 10^ we find 

^l^^Ml = \h.{a,-.Wn{mr^ < 10- X \ln{mr^ . (15) 



^ However, as we will discuss in section II C, this tree level equality can be broken when loop corrections are included. 
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FIG. 2: The exact integral Eq. (12) (solid line) and the approximation Eq. (14) (dashed line) plotted for khoi/kt = \/0.03. 
For this choice, a patch of size £ corresponds to a 3% region of the observable sky, H^^ . On the x-axis wc have plotted 
the values k^^ < Ax < k^^^ — kj^, corresponding to the possible distance between non-overlapping spheres of coordinate 
radius k'^^ /2 within the observable universe. 



For patches £ corresponding to a few percent fraction of the observable universe, the logarithm contributes with 
a factor of order unity, |ln {tH) \ ~ 0(1). ^ The statistical inhomogeneity seen on these scales could therefore be 
at a few percent level. This observation was already made in [46] (see also [47], and the more general discussion 
of [48]). Observational constraints were discussed in [49]. 

It is also interesting to consider the correlators of the quantity ^iC,{-x) = C^(x) — Cjj-i(x), measuring the 
difference in the amplitudes of the curvature perturbation evaluated in the small patch £ and in the full sky 
H~^. To first order in perturbations wc find 

A,C(x) = (7V(x) - (7V(x)),) - (7V(x) - (iV(x))^-i ) (16) 
= (7V(x))^-i - (iV(x)), 

The correlator between A^C(x)'s evaluated on different patches reads 

(AK(xi) A,C(x2))h-i = F{^^^ fchor, h) (17) 

that is proportional to the two point function measured in the full sky. Also, we have 

(AK(xi)A^lnPf(x2))j^-i = -y/NL^^., -F(Aa;,fchor,fc£) (18) 

that is proportional to the three point function. 

In this paper we will focus on the two-point correlators of n-point functions measured in separate patches. It 
would also be possible to extend our analysis to higher point correlators of the same quantities between different 
patches. This would require extending our analysis to non-linear orders in A^cTx and, depending on the model, 
also extending the definition of A^ctx, given in Eq. (4), to non-linear orders in perturbations. This should be 
technically straightforward although it would lead to longer expressions. We will not pursue this here, also since 
the results will be suppressed by higher powers of the power spectrum. Measuring the higher-order correlators 
of observables in different patches would also become computationally more demanding, something which our 
analysis, constraining the higher-point correlators, alleviates. 



B. Correlators of three-point functions 



Because inhomogeneity of the power spectrum is already restricted by observational bounds on the primordial 
bispectrum, it is more interesting to study the correlators of non-Gaussian observables such as /nl itself. The 
non-linearity parameter /nl in a patch (., located at x, can be expressed as 

/NL(a^,x) = /NL(cTff-i) + /^l(<^H-OA^(Tx + 0{f^^\ecjl) . (19) 



* Here and in the rest of this paper, we will always focus on cases in which the size of the logarithms is of order one. The effects 
we are interested in are observable thanks to the coefficients in front of the logarithms, that can assume large values. 
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The prime denotes a derivative with respect to <tjj-i and for single-source models we find 

/NL(^i/-0 = (^SNL(a«-i) - Ja^-i)) N'{aH-^) , (20) 

where ^inl is a parameter describing the primordial trispectrum relative to the power spectrum cubed, and is 
given in the 5A''-approach by [50] 

25 N'" . . 

5NL = ^^- (21) 

For brevity we suppress the argument o-jj-i in what follows, and write /NL(fi/-i) = /nl etc meaning that 
quantities without argument are evaluated in the full sky. 

To leading order in A^CTx, the difference between /nl measured in the entire observable universe and in 
a patch of size £ is then given by 

A^/nl(x) = /nl(x) - /nl = (^^5NL - y/nl) N'Aia^ , (22) 

and the two-point function of A^/nl(x) reads 

'9 12 ,2 ' ' 



(A,/nl(xi)A,/nl(x2)) = (g5NL-y/NLj N"{Aea^,A,a^,) (23) 

9 12 .2 \ ' fHV f''' J , , sin(fc'Aa:) 

S^nl-^/SlJ N-[^) Ij^'^^ 

9 12 \ ^ 

g5NL - y/NL j -F(Ax, /Chor, kg) . 

We then learn that (A^/^l) is proportional to a combination of terms containing powers of fiiNL and /nl- As 
we will discuss in the next sections, those terms are related, by consistency relations valid in the single-source 
limit, to parameters controlling five and six point functions. 

Using the amplitude of spectrum of curvature perturbations, = 2.4 x 10~^, we find 

(A,/nl(xi)A,/nl(x2)) = 7.8 X IQ-'^g^^FiAx, fchor, ke) + O(10-«/^l.9nl) + 0{W-\f^^) . (24) 

For |s(nl| ^ /nl inhomogeneities are unobservably small, ((A^/nl)^) ^ 1, as /nl is constrained by |/nl| ^ 
10^ [9]. On the other hand, if IqnlI ^ /nl' ^® discuss this possibility further in the next subsection, the 
inhomogeneities can become significant as the observational constraint for ^nl is rather weak, |(7nl| ^ 10^ 
[10-14]. In this limit, the last two terms in eq. (24), proportional to powers of /nl, give only subleading 
contributions and they can therefore be neglected. The function F{Ax,khoT,ke) contributes with a factor of 
order unity, and thus we arrive at the result 

(A^/nl(xi)A^/nl(x2)) « 10-«5nl ^ 10^ • (25) 

For IanlI ~ 10^, the variation of /nl measured in different patches of size £ < becomes quite large, 

((A^/nl)^) ^10'*, implying that /nl becomes inhomogeneous on small scales ^. This is a generic feature of all 
single-source models with large ^nl- ^ 

It is illuminating to compare the result with the observational accuracy for detecting /nl using a small 
fraction of the sky. For example, the BOOMERanG experiment found —670 < /nl < 30 at 65% CL observing 
a 3% region of the sky [51]. For patches of the same size, we find the variance ((A^/nl)^) ~ 1-4 x 10~^g^j^ 
using Eq. (24). For I.^nlI ^ 10®, the variation of /nl measured on different 3%-of-the-sky patches is therefore 
of the same order of magnitude as the BOOMERanG 1-a accuracy for detecting /nl- 

In general, when estimating n-point functions on a limited part of the sky, observational errors increase (at 
least) proportionally to 1 / -^//sky (/sky denotes the fraction of the sky which is observed), since there is less data 
available in a small patch. Therefore when dividing the sky into, say, one hundred equal patches, the error bar 
on every observation will grow by at least a factor of ten compared to the full sky measurement. On the other 
hand, we then correlate, over the full sky, the measurements of n and fc-point functions evaluated in distinct 



^ We emphasize that our results are derived retaining only the linear term in (19). For |gNLl ^ /nl' '^hich is the case of interest 
here, the higher order terms 0{fy-^Af(T'^) give a subleading contribution provided that \N /N"^ \ <C IO^IjinlI- 

^ A comparable constraint on the inhomogeneity of /nl can be found by considering the correlator of A^/nl and A^^, see (50) 
for this correlator in the multiple-source case. We thank Antony Lewis for pointing this out. 
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small patches. This implies that, for the aim of constraining the values of higher-order point functions in the 
full sky, we are able to regain much of the accuracy we had previously lost. 

There is however an important second effect to take into account: when making measurements of n-point 
functions in a patch smaller than the entire sky, a smaller number of correlators can be built. For example, for 
n = 3, when estimating /nl in a patch the range of scales available is reduced, thereby limiting our ability to 
analyze the squeezed limit of triangles . This effect is expected to depend only logarithmically on the size of 
the patches. Indeed, the signal to noise ratio for measurements of local non-Gaussianity scales as [52, 53] 

where i'max and ^min denote the largest and smallest multipolcs accessible by the experiment. For example, for 
patches covering one per cents of the sky, fmin/^max -C 1, and the patch size, which determines ^min, enters only 
in the logarithm, reducing the accuracy only by a factor of order three with respect to full sky. We will not 
develop this issue further in our work, but to account for these logarithmic effects we simply increase by one 
order of magnitude the values of our expected constraints on higher order n-point functions in the full sky. We 
will see that, even within this conservative estimate, we can find interesting constraints on these quantities. 

To conclude this section, let us point out that alternative techniques based on needlets analysis of CMB data 
are particularly well suited for testing non-Gaussianity in selected regions of the sky. This fact has already 
been used for investigating inhomogcncity of non-Gaussianity in different, large regions of the sky in [54, 55], 
with the main aim of investigating forcgroimd contaminations and directional variations of /nl • Although no 
significant hints of anisotropics have been reported in those studies, it would be very interesting to apply the 
same techniques to a collection of smaller regions of the sky, to instead test inhomogeneities of /nl along the 
lines we suggest here. 



C. New perspectives for discriminating single from multiple source models 

We have learned from the previous analysis that a large g'NL leads to sizable inhomogeneities of /nl- On 
the other hand, other inflationary observables can be affected by a particularly large ^nl- We discuss in this 
section this possibility, showing that the conditions allowing to obtain large inhomogeneities for /nl offer new 
perspectives for distinguishing single from multiple-source scenarios. 

Loop corrections [56] are known to influence observables by providing logarithmic contributions that in some 
scenarios can be large enough to dominate over tree- level quantities [57, 58]. In our set-up, tnl is the observable 
that is most sensitive to loop corrections when ^nl is large. Expressing it as 

Tnl — Tnl + t^l ' 

one finds that for |(7nl| ^ |/nlP, the dominant part of the one- loop correction to tnl is given by 

gl^V^\x,{^i)<\^\ (27) 



1— loop 
"^NL 



assuming that the fourth order derivative N"" is constrained according to footnote 5. We have used approx- 
imation methods similar to [59 Gl] in evaluating the loop integral. The cut-off scale I should be a bit larger 
than the box in which wc arc making our measurements, such that ln(fc£) > 1: because we arc making a leading 
log approximation one cannot consider the limit of the logarithm going to zero. Comparing equations (23) and 

(27) , we observe that the inhomogeneities of /nl are directly related to the fraction of tnl generated by loops 

((A,/nl)^) ^ tAl'""^ . (28) 

This consistency relation, which holds for I^nlI >• |/nlP, is specific for the single-source case. As we will learn 
in Sec IV A, in multiple source models there are additional quantities contributing to the right hand side of Eq. 

(28) , which in principle make it possible to discriminate between single and multiple source scenarios. 
Although for some models one cannot realise I^nlI ^ /nl' especially for those with quadratic potentials, for a 

review see [62], it is certainly possible to realise |,9nl| ^ /nl in general early universe models. This can happen 
for example in the interacting curvaton scenario [63] , since the value of /nl oscillates depending on the initial 
field value of the curvaton, and the points where /nl = do not coincide with points where 5nl is suppressed. 



We thank Rob Crittenden and Dominic Galliano for discussions on this point. 
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Also an isocurvature field direction during inflation with a quartic potential and zero VEV gives rise to ^nl but 
not /nl [64]. In general finite volume effects will also generate a non-zero /nl in horizon size patches [65, 66], 
but in some patches /nl will still be small while jinl is not. Some models of multifield inflation can also give 
rise to |5nl| ^ /nl ^be end of inflation, see Eq. (82) of [67] (see also [68]). 

Interestingly, large loop corrections to tnl also break the well-known consistency relation tnl = (6/nl/5)^, 
characteristic for single-source models. This relation is indeed not protected against loop corrections, as can 
be nicely demonstrated by considering models where the curvature perturbation takes the form ( = N'S(j) + 
N"'S(j)^ /Q, with no higher order terms. For this class of models /nl = to all orders in perturbation theory, 
within the same patch that this ansatz for ( is valid in. At tree-level tnl also vanishes, t^^^ = 0, but the unique 
one-loop corrections generate a non-vanishing result, tnl ~ (54/25)^.g^L7'^ln(k^) + C'(,9nl^C'' ^ '^^'^^^ 
is no two-loop correction and the three-loop correction is constrained to be much less than unity due to the 
observational bound on ^nl- The relation tnl = (6/nl/5)^ is therefore clearly violated by the loop corrections. 
To the best of our knowledge, this has not been demonstrated previously. However the Suyama-Yamaguchi 
(SY) inequality, tnl ^ (6/nl/5)^, which was originally proved at tree level in [69], and the possible effects of 
loop corrections were discussed in [61, 62, 70], has recently been shown to remain true against loop corrections 
at all orders in a completely model independent sense [71] (although there are subtelties about the definition of 
the local correlators beyond tree level [61]). 

The previous example is by no means a special case: for generic single-source models we also find that the 
relation tnl = (6/nl/5)^ can easily be broken by loops. Assuming there are no large loop corrections with 
higher than three-point vertices, and concentrating on the limit Iqnl] ^ /nl' one-loop corrected /nl is 
given by /nl ~ ./nl''(1 + IS.gNL^c ln(A;^)). Using that [gNL| ^ 10^, we conclude that the loop corrections to 
/nl are small, at most at one per cent level, since IS^nl^c In(fc^) < 10^^. On the other hand, the one-loop 
corrections to tnl of Eq. (27) can be large, significantly altering the tree level relation Tnl'* = (6/nl''/5)^- 
Combining the results so far we find the general one-loop corrected single-source relation between tnl and /nl 
given by 

TNL « (0 ' /Sl + (I) ' 5nl n Hki) . (29) 

The second term on the right hand side arises from the loop corrections and can be as large as 10^, irrespectively 
of the value of /nl- It may therefore give the dominant contribution to tnl- In conclusion, the condition that 
leads to sizeable inhomogcneitics for /nl, that is having a large pnl, also leads to breaking the single-source 
consistency relation between /nl and tnl- On the other hand, even though the loop corrections can cause large 
deviations from the tree-level single-source result t^I'^ = (6/nl''/5)^, the relation between tnl and /nl still 
serves as a useful discriminator between single and miiltiplc soiirce models. Indeed as we will show in section 
IV, the result (27) is converted into a lower bound for Tplj^'""'' in the multiple source case, which also implies 
that (29) becomes a lower bound for tnl- It would be interesting to generalize this analysis allowing for large 
loop corrections with higher than three-point vertices and see if the same conclusion holds even in that case. 



D. Correlators of four-point functions 

Considering correlators involving four-point functions measured in small patches, we can constrain higher- 
order n-point functions in the full sky, with n > 5. In the single-source case, the only new parameter which we 
can independently constrain (even if going beyond 2-point correlators) is the quantity 

1 /5\^7V"" 

This is associated with the fourth order term in the standard expansion of the curvature perturbation 

C = Cg + ^/nlCS+ (0 5nlCS+Q) h^LCh + ---- (31) 



In multiple-source models, more possibilities arise as we are going to discuss in section III; in this subsection, 
we focus on correlators that specifically depend on /inl- These involve measurements of ^nl in small patches, 
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and read ^ 

(A<.C(xi)A|<?nl(x2)) =i 
(A^ln7'(xi)A^5(NL(x2)) 
(A£/nl(xi) A^5(nl(x2)) ^ 

(A^SfNL(xi) AfgNL(x2)) ^ 

where in the right hand side of each correlator we only include the terms proportional to /inl and that are 
weighted by the largest coefficient, under the assumption that the actual values of /nl, and (^nl saturate their 
observational bounds. For simplicity, we also drop the factor F defined in Eq. (12) which is of order unity and 
does not modify our discussion. 

Approximate la constraints with WMAP for logV^, /nl, .9nl are respectively 10"'', 10^^ , 10 , 10®. Con- 
sidering around one hundred patches of the sky, the measurement errors in each patch will be about ten times 
larger than these. On the other hand, when making correlations over the full sky, we should be able to recover 
the accuracy we loose when focussing on a small patch. As discussed at the end of subsection II B, in order to 
take into account effects associated with the fact that only a limited class of momentum space configurations 
fits into the small patches, we increase the expected error bars on the correlations above by one order of magni- 
tude. The resulting error bars are consequently of order 10^, 10®, 10^ and 10^^, respectively for eqs. (32)-(35). 
Considering the constraints these terms could give on /inl, we see the correlation between ( and ^nl, eq. (32), 
would lead to a competitive constraint, of order 

I/inlI <fewx 10". (36) 

Moreover, this correlation is sensitive to the sign of /inl- A comparable constraint can be obtained from the 
autocorrelation of (Jnl, eq- (35), although this quantity cannot determine the sign of /inl ^■ 

Note that the assumption |/inl| IO^I^nlI, which justifies neglecting the second order terms in eq. (19), 
breaks down when I^nl] ~ 10® and the bound |/inl| ^ few x 10^^ is saturated. To work out the precise numerical 
factors in this limit one should therefore include the second-order terms and repeat the analysis leading to eqs. 
(32) - (35). This however should not affect the order of magnitude of the bound |/inl| ^ 10^^. 

How useful is this constraint on /inl? As we did in the previous sections with respect to ^nl, we should 
compare them with constraints obtained from loop contributions to inflationary observables that are sensitive 
to /iNL- There is a "dressed- vertex" loop diagram of the bispectrum [72] which gives ^ ^nl-Pci which 

gives effectively a constraint of the same order as (36) , unless there is a cancellation between loop terms and 
the tree level term contributing to /nl- One would need to include numerical factors to see which constraint is 
the most competitive. 

If one renormalises the 6N coefficients to avoid the dressed vertices, as suggested in [72], one should shift 
the derivative of N introducing N" = N" + N"" {ScjP') /2, which absorbs this term. On the other hand, in this 
way, the fine-tuning needed to compensate the loop corrections does not go away, it just becomes harder to 
spot. If the loops associated with dressed vertices do cancel, then the largest loop term involving /inl becomes 
S'nl'""'' ^ /nl/inl^'c which implies |/inl| ^ 10^^, which is a much weaker constraint than we forecast cound be 
found by considering the autocorrelation of Aj^'nl- 

III. MULTIPLE SOURCE MODELS 

The extension of the previous analysis to the multiple-source case suggests that our method is able to probe 
many of the parameters controlling higher-order point functions. 

Taking ensemble averages of correlators among n-point and fc-point parameters calculated in small patches, 
we are probing parameters associated to {n -{- fc)-point functions evaluated in the entire sky. More precisely, 
for n > fc the correlators of local n- and A;-point functions probe those parameters of the full-sky (n -t- fc)-point 
function that contain at most n:th order derivatives of N ^ the number of e-foldings. This is the generalization 
of the observation made in [71] for the power spectrum. In the single-source limit, some of the {n + fc)-point 
function parameters coincide with the squares of parameters associated with lower point functions (analogously 



Details of these calculations are given by taJiing the single-source limit of the general, multiple-source calculation which we will 
present in section III. 

Planck data will be characterized by error bars reduced by one order of magnitude with respect with the ones discussed above. 
The very same procedure could then be applied, obtaining more stringent limits on /inl- 



^ /nl /iNL 



(32) 
(33) 
(34) 
(35) 
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to the tree-level consistency relation tnl = (6/5/nl) holding between four and three point functions in the 
single source case). This is the reason for which we found that the two-point correlators of A^/nl's calculated 
in different small patches are proportional to the other hand, in the multiple source case, our method 

in principle allows us to probe individually the various parameters controlling the higher-point functions. 

In this section, we discuss this topic in more detail. Recently, a diagrammatic approach, based on the 
^AT-formalism, has been adopted in [16, 73] to determine the parameters characterizing the five and six point 
functions. While the three-point function is characterized by one parameter /nl and the four-point function by 
the two parameters tnl and (7nLj the five and six point functions are characterized by three and six parameters, 
respectively. In this section, we borrow the notation of [16] for defining parameters associated to 5 and 6-point 
functions (see [16] for more details, but notice that we do not adopt their normalization condition NcNc = 1). 



Two point function: 



Three point function: 



Four point function: 



Five point function: 



-Y 



_ 5 NgNbNaj, 



Tnl = 



NgNa^N^bNb 



25 NgNbN.Nabe 

54 {NeNef 



ft 



(1) ^ NqNghNbcNcdNd 



A2) 
J5 



NqNabNbcdNcNd 



f(3) 



NgNhNcNdNgbcd 



(37) 



(38) 



(39) 



(40) 



Six point function: 



(2) 
56 



(NfNff 

NgNgbNbcdNcNdeN^ 
[NfNff 



^(2) ^ NqNbNabcNcdeNdNe 
(3) _ NgNabNbcdeNcNdNe 



(4) 
96 



NgNbNabcNcdNdeNe 

(NfNff 

NgNbNcNdNeNabcde 

(NfNff 



(41) 



(NfNff 

Notice that the parameters associated with 3 and 4 point functions carry the conventional numerical coefficients. 
The aim of this section is to show that ensemble averages of correlators between quantities involving two and 
three-point functions depend on parameters associated with five and six point functions. This suggests a method 
for probing the quantities appearing in eqs. (40) and (41). 

Generalizing the notation of the previous section, we can write in the multiple-field case (summations over 
repeated indices are understood) 

/NL(a^",x) = /NL(a^-i) + /NL,6(a^-i) A^a^ (42) 

where the lower-case upper indices a, 6, ... denote the different scalar fields involved. Analogous expansions hold 
for C, P(^, Tnl and ^nl- We assume that the metric in field space is flat, and that 



i^e<A.<) = (S)' £ 



dk 



, sin (fc' Ax) 
Ax 



(We will neglect the slow-roll suppressed cr-dependence of H throughout this paper). Then we have 
CbAcT^ = -TVfcAa^ 

NahNa , b 



(InP^ , AcT^ 

/nl, d Acr^ 
TNL,d ACTx 
ffNL, d AfTx 



= 2- 



5 {2NadNbNab , NaNbNabd 24 fj.^^ N^N^d 



{NcN.y 



+ 



{NcN.f 



„ / NadNacNcbNb , NaNacdNcbNb „ N^dN^ 
^ I — 1 .ci otnl 



{N,N,f 



{NeN, 



N,N, 



25 I ^ NadNbN^Ngbc ^ NgNbN^Nabcd 324 ^nl iVcrfA^c ] 



54 



{NeNe 



{NeNeT 



25 iVpTVp 



(43) 

(44) 
(45) 

(46) 
(47) 

(48) 
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Denoting 

A^/nl(x) = /nl«,x) - /nl(^t^-i) (49) 

and the same for (, InP^, tnl and ^nl, we obtain the following results (the function F = F {Ax , khoi , ke) is 
given in Eq. (12)): 

(A,C(xi)A,/nl(x2)) = - (^^TNL + ^5NL - y /^l) 'PcF , (50) 

(A,C(xi)A,rNL(x2)) = -2 (^fP + fP ^ ^ tnl/nl) VcF , (51) 

(A,C(xi)A,ffNL(x2)) = - (g /f ^ + I /f ^ - f flNL /nl) , (52) 

(A^lnP(xi)A^lnP(x2)) = W^P^F, (53) 

/A ^ /- ^ A ^ / \\ 25/ (1) (2) , 576/^L , A (1) 96/nl .(i) 48/nl j.(2)\ ^ i-, 

(A^/nl(xi) A^/nl(x2)) = — Ut^ ' +n + TNL + ' — — ' j F ,(54) 

(A,lnP(xi)A,/NL(x2)> = ^ (^4/« + 2/f ) - ^/nltnl) (55) 

(A^ InP(xi) A^rNL(x2)) = 4 (r^^^ + ff^'^ - St^l) ^ . (56) 

(A, InP(xi) A,<7nl(x2)) = ^ f 85^ + 5^ ^ - ^snlTnl) ^ ^ • (57) 



These ensemble averages are able to probe combinations of parameters controlling 5 and 6 point functions, as 
explained above. The previous combinations allow one to individually test each of the parameters characterizing 
the five-point function. Indeed, a combination of Eqs. (51), (52) and (55) lead to the following expressions for 



these parameters evaluated in the full sky 

= I (A,ln7'(xi) A,/nl(x2)> + i(A,C(xi)A,TNL(x2)) + ^/nltnl^c (58) 

/f ) V^F = -/(I) P^F- i(A,C(xi)A,rNL(x2)> + ^/nl tnl P( F , (59) 
(r,\ 1 r9i 54 1944 

/fp^F = --/f)p^i^-_(A,C(xi)A,ffNL(x2)) + ^/NL5NLn^- (60) 



In order to individually test the quantities associated with 6-pt functions one can consider correlations among 
Tnl and ^nl measured in small patches - although the expressions for these quantities depend also on 7 and 
8-pt functions. However the parameter containing fifth order derivatives, can not be probed by considering 

the inhomogeneities of trispectrum to leading order in A^a. In order to probe one should consider the 
inhomogeneities of the 5-point function or include the next-to-leading-order contributions in A^cr. However 
these in general are suppressed by a further factor of "P^. 

It is straightforward to check that, in the single-source limit, one recovers the results of the previous sections. 
This is due to consistency relations associating parameters in eqs (40) and (41) with parameters in (38) and 
(39). In the multiple source case, the parameters of eqs. (40) and (41) are not directly constrained by the 
observational bounds on /nl, tnl and ^nl- The inhomogeneities we are calculating could therefore be very 
large in certain classes of models - even at the level that might be detectable by re-analysing the WMAP data. 

IV. WHEN IS THE INHOMOGENEITY OF NON-GAUSSIAN PARAMETERS LARGE? 

In this section, we investigate in more detail when the autocorrelations of /nl and ^tnl measured in distinct 

patches can be large, and discuss the differences between single and multiple-source scenarios. We focus our 
attention on a representative, two-field case, in which the curvature perturbation is expanded as 

C = + ^Arn502 ^ N^^s^Sx + ^A^22<5x^ + ^A^iii^<^^ + ^Nn2S^^Sx + ^Ni22S^6x^ + ^A^222^X^ -(61) 

An expression of this form can be obtained by doing a rotation in field space to the direction, cf), generating 
the first-order contribution to the power spectrum of (: in this way only one linear term in scalar perturbations 
appear in Eq. (61). One could instead choose a different field basis in which the second order derivatives are 
diagonalised (i.e. A''i2 =0), as was chosen in [70], but this would be at the expense of introducing a second 
linear term which leads to more cumbersome expressions for the non-linearity parameters. We note that, in 
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general, either choice of ehminating one term from the SN expression (61) will lead to a non-trivial field- metric, 
because the required field space rotation will depend on the scale we are focussing in. However since we can only 
observe a limited range of scales (around 5 e-foldings) this should typically only introduce small corrections, 
unless there happen to be large numerical factors which override the expected slow-roll suppression. We will 
not develop this interesting issue further, but simply adopt the expression (61) as a specific Ansatz for the 
curvature perturbations. 

Within this set-up, the non-linearity parameters controlling quantities up to the 6-pt function are: 




A. Inhomogeneous autocorrelation of /nl 

We have learned that, in a single-source set-up, the autocorrelation of /nl can be large provided that we 
saturate the observational limit for the full sky value of ^nl- Let us investigate what happens instead in the 
multiple source case. In order to get an inhomogeneity for the parameter /nl large enough to be potentially 
observable with Planck, we need the combination in brackets of (54) to be of order O{10^'^). 

Given the present observational constraints, |/nl| ^ 10^, tnl ^ 10^, I^/nlI ^ 10^ [9-14], one may check that 

(2) 

the only term which could be so large is the one proportional to Tg , that provides 

(A,/nl(xi)A,/nl(x2)) « |rf Pc^^ = || g^^^) ' + (^^) ' j ^ ^ • (66) 

Notice that the first term in the right hand side with ^fNL is the same leading contribution we found in the single- 
source case, sec Eq. (23). The second term represents a pure multiple-source contribution and its magnitude 
is independent of ^nl- While the form of Tg in Eq. (66) is specific for the two-source model (61) that we are 
considering as an example, the conclusion that the inhomogeneities of /nl are not uniquely determined by the 
magnitude of giNL holds for generic multiple-source scenarios, as we will show below. 

The result (66) shows that large inhomogeneities of /nl are generated by cubic derivatives of N, which also 
generate loop corrections to tnl- In the limit /^^ ^ Isnl], the dominant part of one- loop corrections to tnl is 
given by 



T, 



nZ'°°^ ^{rP + ^ + ^] nin(H) , (67) 



which implies that 



((A./nl)^) - (0 r-'°- - (0 + ^) VMki) < (0 r-'- . (68) 

This has to be compared with what we found in the single-source case, Eq. (28), in which the autocorrelation 
of /nl is uniquely determined by t^l'°°^> equivalently by Q^i^Vq according to the relation (27). In the 
multiple-source case, t^°°'^ instead sets only an upper bound on the inhomogeneities of /nl and the actual 
level of inhomogeneities could be quite different from r^^ 

It is straightforward to show that this conclusion holds for generic multiple-source models, thus opening 
interesting possibilities for discriminating between single and multiple-source scenarios. Defining a unit vector 
wi*^^ = ^ofe, the Cauchy-Schwarz inequality for the inner product of the vectors u^''^ NaNabc and Nc leads to 

{ui'^^NaN^bcNcf < {ui''^NaNabcNcdeNdui''^){NfNf) , (69) 

where only the repeated lower case indices are summed over. The one-loop corrections to tnl containing two 
three-point vertices are given by 

l-loop _ J^a^ubc^hul-^d ^ _ Uf!''' ^a^ubc-^c<U..^dll'!!^ ^ 
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where we have neglected the logarithm associated with loop integrals. This represents the dominant contribution 
to the loop corrections, provided that IqnlI ^ /nl' ^^'^ ^^^^ there are no large loops associated with higher 
than three-point vertices. Similarly we can write 



Comparing the above expressions to the inequality (69), and making use of the inequality (54/25)^ gf^L ^ "^6^^ 
derived in [16], we find 

^^^VcffL<n-f ^riL'""-^. (72) 



.25 

In the limit IqnlI '> /nl' using Eq. (54) we therefore find the general result 



2 

■ 



((A,/nl)^> < ^ r^Z^°°P • (73) 



The inequality is saturated in the single-source case (28) where the magnitude of the inhomogeneities is set 
by t^l'°°'^. In multiple source scenarios, "^^^ greater than ((A^/nl)^)- The inhomogeneities of /nl 

could therefore provide an interesting new tool for discriminating between single and multiple source models. 
While the result (73) only applies in the limit |5nl| ^ /nl assuming there are no large loops with higher 
than three-point vertices, it could nevertheless offer an intriguing new window for probing inflationary physics. 

Similar information can be obtained by considering the structure of the trispectrum. The inequality (72) 
implies 



nS&L < r^L""" , (74) 



which is again saturated for the single-source case. Under the assumptions stated above, there are no significant 
loop corrections to (/nl, and tnl is dominated by the loop correction as we are considering the limit Tnl'' ~ 
/nl ^ \9nl\- We therefore see that both the inhomogeneities of /nl and the ratio of the two parameters 51NL 
and Tnl have the opportunity to provide interesting new tools for distinguishing between single and multiple 
source models. 



B. Inhomogeneities of (;nl 

We can apply the same procedure to analyze the autocorrelator of gfNL- This quantity is given by 

/ X A / ^^ /^625 (1) 625 (3) 2 3750 (2) 50 m 50 (3)\ ^ , 

(Af5iNL(xi) Af5NL(x2)) = ( ^"^8 + ^^^8 " 365nlTnl + ^^98 - y^NL^e - y^NL^e I ^c47.5) 

Based on a measurement accuracy of |(?nlI ~ 10*" (which will be improved by up to two orders of magnitude 
with Planck [14]) we require the term in brackets to be 0(10^^) in order to be able to probe the inhomogeneity 
with WMAP data. With Planck data values as small as C(10^^) could be relevant. 
The new non-linearity parameters which enter at this order are: 

gi^^Nf = NniiNn+Nni2Ni2, Nf = N^,-,, + N^^,^, (76) 

rg^^^ATio = Nfi,Nn{Nn + A^i2) + 2NniNn2Ni2{Nn + N22) + A^n2(^l2 + ^1^), (77) 
5(2) iV9 = 7Vnii(ArniViii+7Vi2A^ii2)+A^iii2(A^i2A^iii+A^22iVii2). (78) 

Notice that only two new quantities are important for studying the new correlator of Eq. (75). These are A^'im 
and A^iii2, which arc two out of the five foiirtli derivatives of N. Unless one of these two quantities is large, the 
correlator (75) does not give competitive information compared to the correlators we considered in the previous 
sections. If one of the fourth derivatives is assumed to be very large, then the dominant term is simply Tg^^ . 
This implies that we can constrain Nnn/Nf and Nni2/Nf to the level of 10^^ with WMAP and 10^^ with 
Planck. 

On the other hand, I/nl'""""! - \Niiii\/NfP^ < 10^ provides the constraint \Niiii\/Nf < 10", so barring 
an accurate cancellation between this term and the tree level /nl° do not find interesting constraints from 
this contribution. However for A^iii2 there is no 1-loop constraint from /nl, and instead the tightest constraint 
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comes from t^l'°°^ ~ -^1112^12^^/^1 ^ This does not constitute a dominant constraint compared to the 
expected sensitivity with the Planck satellite unless \Ni2\/Nf is larger than unity. 

In conclusion, we see some similarities between studying the autocorrelation of A^/nl and here the autocor- 
relation of A^giNL- In both cases one can constrain two higher order derivatives, the one which is present in the 
single-source case, i.e. A''iii and Nun which are proportional to giNL and /inl respectively, but this constraint 
is not very competitive compared to the natural constraint we get from consider the loop contribution to lower 
order correlators unless there is a chance cancellation between the loop and the tree level terms. For the 'cross- 
derivatives', iVii2 and iViii2, the loop constraints are weaker (especially in the case of iVii2) and we can achieve 
a tighter probe by considering the inhomogeneity of non-Gaussian correlators. 

V. CONCLUSIONS 

In this paper, we investigated under which conditions inflationary parameters can be inhomogeneous to an 
observable level, focussing in particular on observables controlling local non-Gaussianity. We have demonstrated 
that if we subdivide the entire sky into a large number of small patches, the value of cosmological observables 
associated with the properties of the curvature fluctuation can differ from patch to patch. In particular, we 
have shown that correlators between n-point functions of curvature fluctuations as measured in one patch, and 
fc-point functions as measured in another patch, depend on (n + fc)— point functions in the entire sky. This 
implies that the expected degree of inhomogeneity in observable quantities can be quantifled in a rather model 
independent manner. In interesting eases it results large enough to be measurable, since at present we have 
only weak constraints on parameters controlling higher-order point functions in the full sky. Consequently, 
inhomogeneities of non-Gaussian parameters can also be seen as feasible method for probing or constraining 
higher-point functions. 

We have analyzed in detail the degree of inhomogeneity of local non-Gaussian observables, flrst in the single- 
source case, in which only one scalar fleld contributes to the generation of primordial curvature perturbation 
(as in curvaton models), then in multiple-source set-ups. In the ease of single-source models, we have shown 
that autocorrelators of /nl evaluated in different patches depend on the value of (?nl in the full sky. If jfNL 
turns out to be large enough to saturate its present day bound, we should expect variations of /nl of order one 
hundred from patch to patch, large enough to be observable. Autocorrelators of .9nl, on the other hand, can be 
used to set constraints on /inL) the parameter characterizing 5-point functions: present day data are accurate 
enough to be able to set an upper bound |/inl| ^ f(3w x 10^^. In the case of multiple source models, we have 
shown that correlators between parameters controlling n-point functions in different patches, with n < 4, can 
be useful for probing individually the several new parameters that characterize five and six point functions. 

We have pointed out interesting connections between the degree of inhomogeneities and loop corrections to 

non-Gaussian observables. Typically, models that lead to sizable inhomogeneities are also characterized by 
large loop corrections to inflationary observables. We have used this fact to determine consistency relations 
between quantities denoting respectively inhomogeneities and loop contributions to non-Gaussian observables. 
These consistency relations take the form of inequalities in multiple source models, that are saturated in the 
single-source limit: consequently, they offer new observational perspectives for distinguishing between single 
and multiple source set-ups. 

The conclusion of our theoretical analysis is that a sizable degree of inhomogeneity in non-Gaussian observ- 
ables is allowed by present day bounds on (;nl and tnl- We have then described observational prospects for 
probing inhomogeneities of non-Gaussian observables, discussing the accuracy we should expect for determining 
correlators of n-point functions measured in different patches of the skies. We discussed geometrical effects 
that reduce the accuracy, and conservatively took them into account in our estimates. We have also pointed 
out that alternative techniques based on needlets analysis of CMB data are particularly well suited for testing 
non-Gaussianity in selected regions of the sky. It would be very interesting to apply those techniques to study 
inhomogeneities . 
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